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Abstract 



o 

I We derive a recurrence relation for the number of simple vertex-labelled 

bipartite graphs with given degrees of the vertices and use this result to 
obtain a new method for computing the growth function of the Artin 
monoid of type A n -i with respect to the simple elements (permutation 
j"y^ braids) as generators. Instead of matrices of size 2 n_1 x 2™-\ we use 

Q matrices of size p(n) x p(n), where p(n) is the number of partitions of n. 

1 Introduction 

Let M be a monoid with a finite generating set S, let S* be the free monoid 
over 5, and let a : S* —¥ M be the natural epimorphism. 

For w £ S* we denote by \w\ the length of w, and for g g M we denote by \g\s 
\jQ \ the length of a shortest word w £ S* satisfying a(w) — g. We call w € S 1 * a 

geodesic representative with respect to S for 5 € Af, if a(w) = g and |w| = \g\s- 
The spherical growth function of M with respect to S is defined as the formal 
power series 
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where is the number of elements of M that have a geodesic representative of 
length £ with respect to S. That is, computing the spherical growth function is 
^ ■ equivalent to counting the elements of given geodesic length. 

It is well known that the spherical growth function is a rational function, if 
M admits a geodesic automatic structure with respect to S |10) ; the latter in 
particular is the case, if each element of M admits a unique normal form which is 
a geodesic representative, and if there exists a finite state automaton recognising 
normal forms. Geodesic automatic structures of this kind have been described 
for the braid monoid [S] and, more generally, Artin monoids of finite type [3] 
and (yet more generally) Garside monoids of spherical type [S] . 

The transition matrix of a finite state automaton with states {1, . . . , N} is an 
N x N matrix T £ Z WxW whose entry Tij is the number of transitions leading 
from state j to state i. If T is the transition matrix of the finite state automaton 
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recognising normal forms, with state 1 being the starting state, u— (1, 0, . . . , 0)*, 
and if v £ Z lxN is the characteristic function of the set of accept states, then 
the number of accepted strings of length £, that is the number of elements of 
geodesic length £, can be computed as ai — vT l u [TU] 

For Artin monoids, two finite state automata recognising normal forms have 
been described in [3]; cf. IScction 3l The states of the first automaton are in 
bijection to the simple elements of the monoid, which in the case of Artin 
monoids correspond to the elements of the associated Coxeter group; this au- 
tomaton readily generalises to Garside monoids of spherical type. The states of 
the second automaton are in bijection to the finishing sets of simple elements; 
the finishing set of a simple element is the set of atoms of M that can occur last 
in any reduced word representing the given simple element. 

In the important case of the Artin monoid of type j4„_i, the braid monoid 
on n strands, these automata have (n — 1)! respectively 2™" 1 states. Due to 
their size, computing the growth function of the braid monoid Z?+ using the 
methods described in [4] becomes infeasible in practice, if n > 10. 

Being able to compute the growth function, respectively the number of elements 
of given geodesic length, is of practical importance, among other things, for the 
purpose of generating uniformly random elements. If an is known, and if the 
elements of geodesic length £ can be enumerated effectively, then a uniformly 
random element of geodesic length £ can be chosen. Such a uniform random 
generator has recently been described for the braid monoid with the set of 
atoms as generating set |11) . (That is, a positive braid of given word length is 
chosen with a uniform probability.) For many purposes, however, being able 
to generate uniformly random braids with given length with respect to the set 
of simple braids would be more useful. This problem is significantly harder; 
improving the methods for computing growth functions is a crucial prerequisite. 

In this paper, we give a new method to compute the growth function of the braid 
monoid B£ that uses a matrix of size p(n) x p(n), where p(n) is the number 
of partitions of n. The basic idea is as follows: We start with the automaton 
from [J] whose states are the finishing sets of simple elements; these finishing 
sets can be identified with the standard parabolic subgroups of the braid group 
B n , respectively those of the associated Coxeter group S„. We will show that it 
is possible to combine states for which the sizes of the orbits of certain subgroups 
of & n agree, arriving at a reduced transition matrix whose rows and columns 
are indexed by the partitions of n. The entries of the latter can be obtained 
from the numbers of simple vertex-labelled bipartite graphs with n edges with 
given degree functions; a vertex of degree k representing an orbit of size k. 
The structure of the paper is as follows: In IScction 21 we derive a recursive 
formula for counting simple vertex-labelled bipartite graphs with n edges. In 
IScction 3l we recall the Artin Garside structure of the braid monoid £?+ , simple 
elements, greedy normal forms, as well as starting sets and finishing sets and the 
finite state automata recognising normal forms described in [3] ; experts may skip 
this section. In IScction 4l we establish that computing the growth function of 
the braid monoid can be reduced to counting simple vertex-labelled bipartite 
graphs with n edges. In IScction 51 finally, we discuss the feasibility of the new 
method. 

After the journal version of this paper was accepted, we learnt that some the 
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results obtained here may be connected to results from [BJ. 



2 Counting vertex-labelled bipartite graphs 

For a vertex-labelled bipartite graph T — (U, V, E), where U — {ity . . . , u r } and 
V = {vi, . . . , v s }, we can consider the sequences degy = [deg(ui), . . . , deg(u r )] 
and degy = [deg(«i), . . . , deg(u s )] which satisfy deg(ui) + • • • + deg(w r ) = \E\ 
and deg(«i) H h deg(u s ) = \E\. 

We want to count the simple graphs yielding two given sequences. For n G N 
and sequences A = [01, . . . , ay] and B = [bx, . . . , b s ] of non- negative integers 
with A = ai H — ■ + a r = n = Y] B = bx + ■ ■ ■ + b s , let Na,b be the number of 
simple vertex-labelled bipartite graphs (U, V, E) with degy = A and degy = B. 

Lemma 2.1. Let A = [cti, . . . ,a r ] and B = [b\, . . . ,b s ], where ^2 A — ^2 B . 

(i) N AfB = N B ,A- 

(ii) If A' and B' are equal to A respectively B up to permutation of entries and 
addition or removal of entries that are equal to 0, then Na',b' — Na.b- 

Proof. Trivial. □ 

Lemma 2.2. If A — [ai, . . . , a r ] and B — [1, . . . , 1], where A = n = B, 
then 

Na^b = — j r ■ 

Proof. The graphs (U, V, E) with deg^ = A and degy = B = [f, .... 1] are in 
bijection to the functions {1, . . . , n} — > {1, . . . , r} for which the image i occurs 
exactly a,j times. Thus, Na,b is the number of permutations of a multiset whose 
sequence of multiplicities is A. □ 

Proposition 2.3. If A — [ay . . . , a r ] and B = [b%, . . . , 6 S ], where Y2,A = Y2,B, 
then 

N A ,B = N AuB- , 

7C{l,...,r} 
1/1=6! 

where B' = [62, ... , b s ] and Aj = [a[, ... , a' r ] with a\ = ai — ( 1/ is the 

indicator function of the set I.) 

Proof. Consider the vertex v\ of degree b\. As we want to count simple graphs, 
there can be at most one edge from V\ to every vertex in U; the b\ edges ending 
at V\ thus determine a unique oysubset of {1, . . . , r}, and vice versa. 
For a fixed configuration of edges ending at v\ , the remaining n — b\ edges form 
a simple vertex-labelled bipartite graphs (V, V', E'), for which degy/ = Ai and 
degy, = B' with Aj and B' as in the statement of the proposition. Thus the 
claim is shown. □ 

Bv lLcmma 2.11 Na,b only depends on the types of the partitions of n given by 
A respectively B, and we have Na.b — Nb,a- There are thus hp(n)(p(ri) + l) 
numbers to compute, where p{n) is the number of partitions of n. 
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Corollary 2.4. The numbers Na,b for ^2 A — ^2 B < n can be computed in 
time O(c^) with c < 569. 

Proof. First assume that Na'.b' is known and bounded by f(n — 1) for all 
sequences A' and B' with ^ A' — ^ B' < n — 1. 

For S — [si, . . . , St] define min(5) = min {si,..., s t } and \S\ =t. By[L cmma 2Tl 
we may assume that b\ = min(-B) < min(A) in the statement of |Proposition 2?3j 
and thus \A\ < The equality for Na.b therefore involves 

^<(^A<( n 4\ bl <e> 



h J ~ V h J - \ b\ 

terms, whence Na,b can be computed in time 0^e 2 V^ ln/(n — 1)^. Moreover, 
In /(re) < In/(n-l)+2 v /f . As /(l) = 1 bv lLemma 23\ one has ln/(n) < n 2 . 
Thus, knowing Na'.b' for all sequences A' and £?' with ~^2 A' ='^2 B' < n — 1, 
the numbers JV^s for all sequences A and B with ^ A = J2 £> = n can be 

computed in time o{\^\ 2 n 2 e 2 ^^^j = 0(c^™) with c < 569, where the last step 
follows with the asymptotic expression [TJ] 



As 2m=i < ricV™' and n € O(a^) for any a > 1, the claim follows. □ 

Remark 2.5. The choice 6i = min(£>) < min(yl) which was used in the proof of 
|Corollary 2.4| need not yield the smallest possible number of terms. In practice, 
one should compare the binomial coefficients 

( lAl ) ( lAl ) ( ]B] ) and ( |B| 
\min(_B) J ' \max(£>) / ' \min(^4) J \ma,x(A) 

where max ([si, . . . , s*]) = max{si, . . . , s t }, and apply [Proposition 2.3| to the 
situation that corresponds to the smallest binomial coefficient. 



3 Garside structure and geodesic automation 

In this section, we recall the finite state automata recognising normal forms that 
were described in [4], using the terminology of Garside monoids [5 . For details, 
we refer to [TJ [31 HI [SI [3 [HI [H] - Experts may skip this section. 

3.1 Garside normal forms 

Garside monoids of spherical type were shown to be automatic in [5]. The 
automatic structure of a general Garside monoid of spherical type is closely 
related to the greedy normal form; it is a generalisation of one of the automatic 
structures described in [4]. 

In a cancellative monoid M, we can define the prefix partial order: For x, y G M, 
we say x =^ y if there exists c € M such that xc = y. Similarly, we define the 
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suffix partial order by saying that x I== y if there exists c £ M such that x = cy. 
We call seMan atom, if s = ab (with a, b 6 M) implies a = 1 or 6 = 1. 

A cancellative monoid M is called a Garside monoid of spherical type, if it is a 
lattice (i.e., least common multiples and greatest common divisors exist and are 
unique) with respect to =^ and with respect to >p, if there are no strict infinite 
descending chains with respect to either =<; or )>=, and if there exists an element 
A £ M, such that D = {s£M\s4A} = {s£M\A^ s} is finite and 
generates M. In this case, we call A a Garside element and the elements of D 
the simple elements (with respect to A). Moreover, we denote the ==!-gcd and 
the ^-lcm of x, y £ M by x A y respectively x V y, and the ^=-gcd and the ^=-lcm 
of x, y £ M by x A y respectively x V y. It follows that, for s £ D, there exists a 
unique clement ds £ D such that s ds = A. 

We assume for the rest of this section that M is a Garside monoid as above. As 
D generates M , every element x £ M can be written in the form x = X\ ■ ■ ■ Xk 
with x\, . . . ,Xk £ D. The representation as a product of this form can be made 
unique by requiring that each simple factor is non-trivial and maximal with 
respect to =<;. More precisely, we say that x = x\ ■ ■ ■ x% is in (left) normal form, 
if Xk ^ 1 and if xi = A A(xi ■ ■ ■ Xk) for i = 1, . . . , k. Equivalently, we can require 

(1) Xk =/= 1 and dxi A a^+i = 1 for i = 1, . . . ,k — 1. 

Since M is cancellative, it follows easily by induction that the normal form of 
x G M is a geodesic representative of x with respect to the generating set D. 
Importantly, (1) is a local condition, only involving pairs of adjacent factors. 
Moreover, if X\ ■ ■ ■ Xk £ D* is in normal form, then so is every initial subword 
x\ ■ ■ ■ X£ for £ < k. Thus, we can check whether a word x\ ■ ■ ■ x% £ D* is in 
normal form by reading it from the left to the right, only keeping track of the 
last read symbol: If x\ ■ ■ ■ x^ is in normal form, testing whether xg+\ ^ 1 and 
dx£ A X£ + i = 1 suffices to decide whether x\ ■ ■ ■ X£ + i is in normal form. 

In other words, the following finite state automaton recognises normal forms: 



states 


initial state 


accept states 


transition function 


D 


A 


D\{1] 


\xi \idX/\Xi = l 
1 \l if dX A x t £ 1 



Note that, since dA — 1, the automaton will be in an accept state after reading 
the first symbol x\, unless x\ — 1. (A single non-trivial simple element is in 
normal form.) Note further that, since dl — A, once the automaton has reached 
the fail state 1, it will remain in this state. 



The Garside structure of Artin monoids was described in [31 [7]. If M is an 
Artin monoid, its simple elements are lifts of the elements of the corresponding 
Coxeter group W , its Garside element is the lift of the longest element of W , 
and its atoms are lifts of the simple reflections in Moreover, the Garside 

^^We remark that Artin groups, apart from the Garside structure used here, admit another 
Garside structure, the so-called dual Garside structure [T; the Garside structure used here 
is known as the classical Garside structure. However, the monoids (of positive elements) 
associated to these two Garside structures of an Artin group are different. The monoid we 
refer to as "Artin monoid of type T" is the monoid of positive elements of the Artin group of 
type T with respect to the Coxeter- type generators. 
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element is equal to both the =^-lcm of the atoms and the ^=-lcm of the atoms. 
In the case of the Artin monoid of type A n _i, the braid monoid on n strands, 
the set of simple elements is in bijection to the symmetric group ©„_i, so the 
automaton arising from the Garside structure has (n — 1)! states. 

3.2 Starting and finishing sets 

In the special case of Artin monoids, an element x G M is simple if and only 
if it is square-free, that is, if and only if it cannot be written as x — ua 2 v with 
a,u,v G M and a ^ 1. This property makes it possible to recognise normal 
forms with an automaton involving fewer states. We assume for the rest of this 
section that M is an Artin monoid with Garside element A and set of simple 
elements D. Let A denote the set of atoms of M. 

Given x G M, we define the starting set of x as S(x) = {a G A \ a =4 %} and 

the finishing set of x as F(x) — {a G A \ x ^= a}. For any x G D, we have 

F(x) n S(dx) = and F(x) U S(dx) = A 0J Lemma 4.2]. For u,v G D, one 

therefore has du A v = 1 if and only if S(v) C .4. \ S(du) = F(u). 

We can thus check whether a word a;i • • • G D* is in normal form by reading 

it from the left to the right, only keeping track of the finishing set of the last 

read symbol: If x% ■ ■ ■ X£ is in normal form, testing whether x^+i ^ 1 and 

S(xe+i) C F{x() suffices to decide whether x\ ■ ■ ■ x^+i is in normal form. 

In other words, the following finite state automaton recognises normal forms; 

^(A) denotes the power set of the set A of atoms: 



states initial state accept states transition function 



Note that F(xi) = if and only if Xi — 1. Thus, the automaton will be in an 
accept state after reading the first symbol x\, unless X\ = 1. Moreover, once 
the automaton has reached the fail state 0, it will remain in this state. 

In the case of the Artin monoid of type A„_i, the braid monoid on n strands, 
the set of atoms has n — 1 elements, so the automaton described in this section 
has 2" _1 states. 

Example 3.1. Consider the braid monoid 

B 4 = ( 0-1,(72, (73 | (71(7201 = O2O1O2, = O^C^Cs, ^1^3 = 03(7i) . 

The set of atoms is A = {171,(72,03}, and the simple elements are in bijection 
to 64. Thus, the automaton has 8 states, and there are 24 transitions originating 
from any given state. The states and the transition matrix are given below. 

Observe that the rank of T is 5. Columns 1, 3 and 8 of T are unique, but the 
following columns of T are identical: columns 2 and 5; columns 4, 6, and 7. 



V(A) A 
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state finishing set 



1 





2 


Wi} 


3 


W 


4 


{0-1,0-2} 


5 


W 


6 


{0-1,0-3} 


7 


{0-2,03} 


8 


A 



T : 



/24 


21 


19 


13 


21 


13 


13 


A 





1 


1 


2 


1 


2 


2 


3 





1 


2 


3 


1 


3 


3 


5 











1 





1 


1 


3 





1 


1 


2 


1 


2 


2 


3 








1 


2 





2 


2 


5 











1 





1 


1 


3 


V 




















V 



4 Growth function of the braid monoid 

For an Artin monoid with set of atoms A we define \/X = xi\/---\/xkED 
for X = {xt, . . . , x k } C A and V = 1 G D. We analogously define \/ ^- 

Lemma 4.1. Let M be an Artin monoid with set of atoms A and set of simple 
elements D, and let X, Y C A. Then 

{s G D I S(s) C X and Y C F(s)} 

|t6D I t(V^) A (\M\^) = 1 a « rf *A(Vy) = l} 

Proof. The submonoid of M generated by Y is an Artin monoid with Garside 
element El- I n particular, \f Y = \J Y . As an element of M is simple if 
and only it is square-free, this implies that t( \/ Y) , ior t £ D, is simple if and 
only if no atom in Y is a suffix of t, that is, if and only if t A ( V = 1. 
For s G D one has S(s) C X if and only if s A (V -A \ X) = 1, and Y C F(s) if 
and only if s ^ \j Y . Moreover, s )^\JY with s G I? is equivalent to s = t( Y Y") 
for t € D such that t( V Y") is simple, that is, such that tA[\/Y)=l □ 

For the remainder of the paper, we restrict our attention to the braid monoid + . 
The set of atoms is A = {ci, . . . , o n _i} and the canonical map n : D — > & n is 
a bijection, where 7r(er,) is the transposition (i i + 1). In particular, s G D is 
uniquely determined by the permutation 7r s = 7r(s) € © n it induces on the n 
strands. More precisely, for s *E D and i — 1, . . . , n — 1, one has <Xj ^1 s if and 
only if ir s (i) > n s (i + 1), and s crj if and only if 7r7 1 (i) > 7r7 1 (i + 1) [9]- 

For any S A, the subgroup (tt(S)) of ©„ generated by 7r(5) is a standard 
parabolic subgroup of 6„. Each orbit of (7r(5)) is a set of consecutive points; 
an orbit of (tt(S)} of size s > 1 corresponds to a maximal irreducible standard 
parabolic subgroup of (tt(S)) generated by s — 1 transpositions. Let Orb(S') = 
[Oi, . . . ,O r ] be the orbits of (n(S)), arranged in their natural order (that is, 
such that for 1 < a < b < r, one has i < j for any i G O a and any j G Ob), and 
let deg(S) = [\Oi\, \O r \] be their sizes. 
The above immediately yields the following lemma. 

Lemma 4.2. For any sequence A = [a\, . . . , a r ] with ^ A = n, there exists a 
unique subset S G V(A) such that deg(S') = A. 
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Proof. The sequence deg(S), together with the fact that the orbits are arranged 
in their natural order, determines the orbits of (tt(S)}. Moreover, the points i 
and i + 1 are in the same orbit of (^(S 1 )) if an d only if Ui € S. □ 

Recall from lScction 21 that for two sequences A and B of non-negative integers, 
each of whose sum is n, we denote by Na,b the number of simple vertex-labelled 
bipartite graphs (U, V, E) for which degy = A and degy = B. Further, recall 
the automaton recognising normal forms described in IScction 3.21 we call this 
automaton X. The states of % are subsets of A, namely the finishing set of the 
last read symbol Xi € D. 

Theorem 4.3. If A, Y C A and Y ^ 0, then the number of transitions from 
the state Xofl-to an accept state Y' of % with Y C.Y' is N^ eg (y) t deg(A\x)- 

Proof. A transition s from the state X ends in an accept state if and only if 
S(s) C X. Bv lLemma 4.11 it only remains to show that 



Ndeg(Y),deg(A\X) 



\teD | t(VV) A (\JA\X) = 1 and t A (\JY) = l\ 



We consider vertex-labelled graphs with vertex set Orb(y^l \ X) U Orb(Y). Let 
Ovb{A\X) = [Oi,...,O r ] and Orb(Y) = [P U ...,P S ]. Given t e D with 
t(\JY) A (\JA \ X) = 1 and t7\(\/Y) = 1, we define a graph T t with vertex 
set Orb(„4 \ X) U Orb(Y) as follows: For i = 1, . . . , n we add an edge from O a 
to P;,, where O a is the element in Orb(.A\ A) containing i and Pb is the element 
in Orb(Y) containing 7r t (i); it is obvious that this yields a bipartite graph and 
that the vertices have degrees deg(.A\ X) respectively deg(Y). 

Let u = \JY and s = tu. As t A (\JA \ X) ^ s A (\JA \X) =1, we have 

i"t(i) < 7i"t(j + 1) if {j, 3 + 1} C O a for some a. Similarly, as t A (V^) = 1) wc 
have 7rj (fe) < TT^~ 1 (k + 1) if {/c, + 1} C A for some &. If there were i, i' € O a 
with i < i' and 7r* (i), TTt(i') € Pfe for some a and &, the above would imply 
+ 1, . . . ,«'} C O a and {% t (i) , <K t (i + 1), . . .,ir t (i')} Q Pb, with ir t (i + 1) = 
7Ti(i) + 1. In particular, aw^ =<; w, and then 7r s (i + 1) = 7r u (7Tt(« + 1)) = 
7Tu(7rt(«) + 1) < 7r u (7r t (i)) = 7r s (i), contradicting s A (\/A \ X) — 1. We have 
thus shown that the graph T t is simple. 

Conversely, let T be a simple vertex-labelled bipartite graph with n edges and 
vertex classes Orb(*A\A) and Orb(Y). We construct 7r tr g ©„, and thus tr € £>, 
by defining 7Tt r (i) for i = 1, . . . ,n (in this order): If i € O a , then let 7Tt r (i) be 
the smallest number in the union of the vertices adjacent to O a that has not 
been chosen as an image. If + 1} C O a for some a then by construction, 
n r (j) < n trU + !)) implying t r A (V-4. \ X) = 1. Similarly, if {k, k + 1} C P b 
for some 6 then by construction 7r^, 1 (/c) < vr f T 1 (fc + 1), implying f A (V^) = 1- 
Finally, let u = \/Y and s = tru. The fact that T is simple implies that, if 
{ji3 + 1} Q O a for some a, then 7Tr t (i) an d TTTtO' + 1) > ^r t (j) lie in different 
orbits of (tt(Y)), whence 7r s (j' + l) = 7r u (7Tr t (j + l)) > n u (irr t (j)) = ^(j)- Tnus > 
s A (V-4\ A) = 1. 

It is obvious from the construction that one has Ft r = T for any simple vertex- 
labelled bipartite graph T with n edges and vertex classes Orb(.4 \ A) and 
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Orb(Y). Conversely, let t u t 2 G D such that U(\/Y) A (\JA \ X) = 1 and 

< S :A (\/y) = 1 for i = 1,2, and assume T tl = T t2 . The latter implies that 
7Ttj = 7r t2 up to permutations of the points within each orbit (or vertex). More 
precisely, there exist a € (7r(.A \ X)) and /3 e (tt(Y)) such that 7r t2 = air^P- 
However, if O a = {j, j + 1, . . . , j + k} and P& , . . . , Pb fc are the vertices adjacent 
to O a , where bo < ■ ■ ■ < bk, then ti A (\/A \ X} = 1 for i = 1,2 implies 
that 7rt 4 (j + £) € Pb f for I = 0, . . .,k and i = 1,2; that is, a = 1. Similarly, 
£j A (\/y) = 1 for i = 1,2 implies /3 = 1, that is, 7r tl = 7r t2 and thus ti = t 2 - We 
have then shown that the maps t H >• r t and r i— » ir; are inverses of each other, 
completing the proof. □ 

Let T be the transition matrix of T; for X,Y Q A the entry Tyx 1S the number 
of transitions from the state X of % to the state Y of X. 

Corollary 4.4. If Xi, X 2 C „4 smc/i i/iat </ie partitions ofn given by deg(A\Xi) 
respectively deg(^4 \ X2) are egita/, i/ien Tz,x x — Tz,x 2 f or a ^ Z C A, that is, 
the columns X\ and X 2 of the transition matrix of % are identical. 

Proof. BvlLcmma 2.ll we have iVdcg(y),dcg(.A\x 1 ) = ^dog(y) : dc g (^\x 2 ) for all 
Y C A. Applying ITheorem 4.31 for all 7^ Y C A and using Mobius inver- 
sion [13] for the poset (V(A), 2), one obtains Tz t Xi = Tz,x 2 f° r all Z 7^ 0. As 
the sum of each column of T is n!, the claim for Z — follows. □ 

Remark 4.5. It is not true in general that Tz 1 .x = Tz 2 .x, if the partitions of n 
given by deg(Zi) respectively deg(Z2) are equal. The reason is that the equality 
of partitions does not necessarily extend to all supersets of Z\ respectively Z 2 . 

As an example, consider Z\ = {cti} and Z 2 = {o~ 2 } for Bf with A — {01, o~ 2 , (T3}. 
One has deg(Zi) = [2, 1, 1] and deg(Z2) = [1, 2, 1], so the partitions of 4 that 
correspond to deg(Zi) and deg(Z 2 ) are equal. However, for Z' — {01,(73} 2 
Z\ one has deg(Z') = [2,2], but no superset of Z 2 yields this partition of 4: 
deg({(7i,o 2 }) = [3,1], deg({a 2l a 3 }) = [1,3], and deg(A) = [4]. 

That is, the symmetry of A r c j eg (y),deg(^\x) under resorting of deg(Y) is broken 
by the Mobius inversion. 

Example 4.6. Consider again the braid monoid Bf . The following table lists, 
for each state, the orbits of the subgroup of 64 generated by the complement 
of the finishing set, and the partition of 4 given by the sizes of these orbits. 



state 


F 


A\F 


Orb(^ \ F) 


partition of 4 


1 





A 


[{1,2,3,4}] 


[4] 


2 


{01} 


{02,03} 


[{1};{2,3,4}] 


[3,1] 


3 


{02} 


{oi,cr 3 } 


[{1,2}; {3, 4}] 


[2,2] 


4 


{0-1,0-2} 


{03} 


[{i};{2};{3,4}] 


[2,1,1] 


5 


{03} 


{0-1,0-2} 


[{1,2, 3}; {4}] 


[3,1] 


6 


{01,03} 


{02} 


[{1};{2,3};{4}] 


[2,1,1] 


7 


{02,03} 


{01} 


[{1,2}; {3}; {4}] 


[2,1,1] 


8 


A 





[{1};{2};{3};{4}] 


[1,1,1,1] 
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Observe that the columns of T (cf. |Examplc 3.1[ ) corresponding to two states 
are identical, if (and only if) the partitions of 4 given by the complements of the 
finishing sets of the states are equal. Note also that there are no three identical 
rows of T, although, up to reordering, deg(Xi) — [2, 1, 1] for each of the three 
sets Xi = {di} for i = 1, 2,3 (cf. IRcmark 4.5[) . 

The idea now is to group states X G V(A) into classes corresponding to the 
partition of n given by deg(^4 \ X), and to work with a matrix whose rows and 
columns are indexed by these classes. 

Notation 4.7. 

(i) Let ^3 denote the set of partitions of n. Let ip : V(A) — > *P be the map 
assigning to each X G V(A) the partition of n given by deg(_4 \ X), and 
choose ip : CP — >■ P(A) such that ip(ifi(a)) = a for all a G 

(ii) Define T G I? 5 *^ by setting Tp >a = ?V^(q)- 



By \Corollary 4-4 Tp iCt does not depend on the choice oftp. 



(hi) Define P G Z^ xV( - A ^ and Q G Z v( - A ^ x< $ as follows: 

P a .x = \\ ij f< ? respectively Q x , a = < * ^ X 

otherwise (J otherwise 



(iv) Define N G lP x<9 as follows: 



N, 



n\ i 



a, p 



if a =[!,...,!] 



N a ,p otherwise 
(v) Define C G Z v ^ xV ^ and S G ^(A^riA) as f Mows: 

fl ifX=A\Y fl ifXCY 

Cx,y = < n ,, . respectively S x ,y = < „ 

I U otherwise I (J otherwise 

Proposition 4.8. 

(i) The matrix S is invertible and S" _1 is the matrix given by 



(S 



)X,Y = 



(_i)in*i ifxcr 



othe 



(ii) One has T = S~ 1 CP t NP and T = PTQ = PS^CP'N. 
(hi) Let uefW xl and v€Z lxV ^ be defined by 

fl ifX = A , [O ifX = % 

ux — \ respectively vx = \ 

otherwise 1 otherwise 
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and let u € Z^ xl and v € Z lx< ^ &e defined by 
fl i/a = ^) = [!,...,!] 



otherwise 



and v 



i/a = <p(0) = [n] 

1 otherwise 



Then for any k £ N, one has vT k u — vT k u. In particular, the growth 
function of can be computed using the matrix T. 

Proof. (i) The claim is just a restatement of Mobius inversion [13] for the 
poset (V(A),2)- 

(ii) It is elementary to check that T = PTQ, that PQ = 1 e Z^ x<v , and that 



{CP t NP) Y ,x = 



N, 



dcg(Y),dcg(^\X) 



if Y = 
otherwise 



Bv lTheorerriT3l and part (i), the latter implies T = S^ 1 CP t NP proving 
the claim. 



(hi) Since ^([n]) = {0} and 
tively Qu — u. Since 



(QP) 



X.Y 



, 1]) = {-4}, one has vP = v respec- 

1 i£X = 1>((p(Y)) 
otherwise 



and Tx^( v (y)) = T x ,y for all X,Y € V(A) by |Corollary 4.4[ one has 
TQP = T. Thus, for any k e N, 



uT fc 2 = v{PTQ) k u = (vP)T k (Qu) = vT k u . 



□ 



Example 4.9. Consider again the braid monoid BJ ; recall Example 3.1 and 
Example 4.6) With = {[4], [3, 1], [2, 2], [2, 1, 1], [1, 1, 1, 1]}, one has 



N 



P = 



fo 













1\ 

















1 


4 














1 


2 


6 




1 





1 




2 


5 


12 






V 


4 




6 


12 


24/ 






fi 




















°\ 





1 








1 

















1 


























1 





1 


1 

























V 



N 



( 











1\ 











1 


4 








1 


2 


6 





1 


2 


5 


12 


\24 


24 


24 


24 


24/ 




fl 








°\ 







1 
















1 













1 





Q = 












































Vo 








V 
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c = 



/n 


A 

u 


U 


A 

u 


A 

u 
















/1 1 


1 


1 1 


1 


1 

1 


1 \ 


A 
U 


A 

u 


a 

u 


A 

u 


A 

u 





1 


A 

u 








U 1 


U 


1 


-1 

1 


A 

u 


1 

1 


A 

u 


u 


a 
U 


A 

u 


u 


1 





A 

u 








a a 
U U 


1 


1 


A 

u 


1 

1 


-1 

1 


A 

u 


u 


U 


A 

u 


1 

1 








A 

u 






A A 
U U 


A 
(J 


1 


A 

u 


A 

u 


1 


A 

u 


a 
U 


a 
U 


1 

1 


A 

u 








a 

u 




s = 


A A 
U U 


A 
(J 


1 


-1 

1 


1 

1 


1 


A 

u 


a 
U 


1 


u 


u 








A 

u 








A A 
U U 


A 
(J 





-1 

1 


A 

u 


1 





1 






































1 


1 


V 1 




























\0 














V 










/-I 


-1 


-1 




1 


-1 


I 


1 

1 


-1\ 





















1 









-1 





— 1 





1 


















(J 





1 




-1 








—l 


1 














s- 1 




o 










1 





o 


o 


-1 
































1 


-1 


-1 


1 



































1 





-1 






































1 


-1 








































V 






















/24 


21 


19 


13 


1\ 































2 


2 


4 


6 






















f 









1 


2 


3 


5 


































1 


4 


11 


























V 














V 















Our objective is to compute T directly, without summing explicitly over the 
elements of V(A). In the light of proposition 4.8| (ii), the remaining problem is 
to compute PS~ 1 CP t £ Z^ x ^ efficiently that is, to perform Mobius inversion 
using classes of states rather than individual states. 

Notation 4.10. We say that X £ ViA) is of partition type t(X) — a, if the 
partition of n given by deg(X) is a. (That is, t(X) = <p{A \ X).) 
For a £ let N a be the number of subsets X £ ^(A) of partition type a. 
Moreover, let M £ Z*P X ^ ; where M a ^ is the number of subsets Y C X of any 
fixed subset X of partition type a, such that t(Y) — f3. (It is obvious that M a ^ 
does not depend on the choice of X.) 

Example 4.11. For n = 4, consider a — [2,2] and (3 = [2,1,1]. One has 
M Qij 3 = 2, since X = {171,(73} is of partition type a, and X has exactly two 
subsets of partition type /3, namely {o~i} and {173}. 

Proposition 4.12. Let -{{ai*mi, . . . , a r *m r ]t denote a multiset with pairwise 
distinct elements a±, . . . ,a r whose multiplicities are iiifli) = mi for i = 1, . . . , r. 

(i) If a £ and a — §ai*mi, . . . , a r *m r ]j- as a multiset, then 



mi! • • • m r \ 

(ii) For any m £ N, let *}5 m denote the set of partitions of m and define 
Q m = §7*iV 7 I 7 £ CPmJ- with N-y as in (i). Moreover, form\, . . . , m s £ N, 

let K mi ,....m s '■ Qmi x ■ ■ ■ x Q ms — > mi -| hin 3 ^ e ^ e niap induced by the 

joining of partitions. 
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If a = [a\, . . . , a r ] G then M a ^ is the multiplicity of ft in the multiset 

K a 1 ,...,a r {Q-a 1 X • • • X Q ar ). 

(iii) The matrix M can be computed in time 0(d v '™) with d < 2198. 

Proof. Claims (i) and (ii) are obvious from ILcmma 4.2l For claim (iii) consider 
a = [ai,...,a r ] and compute M a ^ for all (3 G ^3 as follows. Start with the 
multiset = -{{ [ ]}}•, where [ ] is the partition of 0; then for i = 1, . . . , r compute 
Q' = Ka.ai (flx£2 a .) (where a = ai+- • - + ai_i) and replace by £2'. Computing 
A 7 for all 7 G <p ai has a total cost of at most 0(\^\n 2 (\nn) 2 ) C 0(|q3| 2 ). At 
any time, £3 and Q ai each contain at most different elements, so there are 
at most |*$| 2 joins to compute in each step, each at a cost of at most 0(n). All 
multiplicities are bounded by 2™, so the product or the sum of two multiplicities 
can be computed at a cost of at most 0(n 2 ). Thus, the cost of each step is at 
most 0(|*P| 2 (1 + n + n 2 )) = 0(n 2 |*p| 2 ). The number r of steps is bounded by n. 
Since there are |CP| possibilities for a, the matrix M can be computed in time 
0(n 3 |*p| 3 ). The claim then follows using the asymptotic expression [T2] 

p{n) 



Proposition 4.13. For a = [oi, . . . , a r ] € <£ ; let sgn(a) = (-l)"" r . Then 
(PS^CP*)^ = sgn(a) sgn(/3) £ sgn( 7 )A 7 M 7 , Q M 7i/3 

for a,/9 G <P- In particular, (PS~ 1 CP t ) a = {PS^CP 1 )^ q . 

Proof. Recall that for A G 'P(.A), the orbits of (7r(A)) correspond to the 
maximal irreducible standard parabolic subgroups of (ir(X)). More precisely, 
an orbit of size s > 1 is generated by the image of a maximal set of con- 
secutive atoms in A that has size s — 1. Hence, if deg(A) = [a\, . . . ,a r ], 
then \X\ = Y.\=i( a i - I) = n - r. In particular, sgn(r(A)) = (-l) |x| = 
(-1)1-^1 (-1)1^1 = sgn(p(X)). 

Let a, (3 £ be fixed. For GeV{A) define 



n G = 

and 

n G = 



{ (A, Y) eV(A)x V{A) | t(A) = a, r(F) = (3, X U Y = G) 



{ (A, Y~) G V(A) x -P(yl) I r(A) = a, t(Y) = (3, X U Y C G} 



If t(G) = 7 then = M 7iQ ,M 7j ^. Moreover, by Mobius inversion [13] for the 
poset (V(A), C), one has 

n G = (-1) |GXG '' "G- ■ 

G'e-P(^) 

G'CG 
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With proposition 4.8| (i) and the above, we obtain 



Xe^ 1 (a) X'eT~ 1 (a) 

Ye v -\p) Y'er-\f)) 

XCA\Y X'UY'=A 



= (-!)" sgn(a) sgn(p)n A 



= (-1)" sgn(a) sgn(/3) ]T (-1)^ G 'I n G , 

= sgn(a) sgn(^) Y N i s S n (l) M 7iCe M 7i(3 

as claimed. □ 

Corollary 4.14. The matrix PS" 1 CP 1 can be computed in time 0(d^) with 
d < 2198. 



Proof. The claim follows from |Proposition 4.13| and |Proposition 4.12"| (iii), noting 
that all entries of M are bounded by 2™, and that n 2 |Cp| 3 e 0(d^ n ) . □ 



Example 4.15. Consider the braid monoid Bf; cf. |Examplc 3.1[ Example 4.6 
and |ExampleT9l With «P = {[4], [3, 1], [2, 2], [2, 1, 1], [1, 1, 1, 1]}, one has 



M = 



(I 


2 


1 


3 


1\ 




(-1 


2 


1 


-3 


1\ 





1 





2 


1 




2 


-2 


-2 


2 











1 


2 


1 


PS" 1 CP 1 = 


1 


-2 





1 














1 


1 




-3 


2 


1 








\o 











V 




V 1 











0/ 



5 Feasibility in practice 

While its time complexity and its space complexity are still exponential in n, 
the new method presented here allows to compute the growth function of braid 
monoids for much larger values of n than the methods from [3]. 

ITable ll shows the time requirements for both methods. Computations were done 
with a development version of Magma [2] V2.18 on a GNU / Linux system with 
an Intel E8400 64-bit CPU (core: 3GHz, FSB: 1333 MHz) and a main memory 
bandwidth of 4.7GB/s (X38 chipset, dual channel DDR2 RAM, memory bus: 
1066 MHz). The maximum amount of memory used was 1370MB. 

Computing the transition matrix of the more efficient of the automata from [4] 
(cf. ISection 3.2p involves considering n! transitions for each of the 2™" 1 states, 
so has time complexity 0(n\ 2™); the space required is 0(4"). On our computer, 
the computation time would well exceed 10 hours for n = 11; the transition 
matrix would exceed 4 GB for n = 16 and 2EB (2 • 10 18 B) for n = 30. 
The new method described in this paper, on the other hand, has time complexity 
O(d^) C 0(2") and requires space 0(|*P| 2 ) = 0(n~ 2 f^) C 0(2") with 
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n 


7 


8 


9 


10 


15 


20 


25 


28 


29 


30 




0.6s 


8.5s 


145s 


2820s 
















<0.1s 


<0.1s 


<0.1s 


<0.1s 


<0.1s 


0.53s 


17.1s 


95.8s 


191s 


359s 



Table 1: Time required for the method from [3|, and time i[ G ] required for 
the new method. 

/ < 170. While the methods from [3] become infeasible if n > 10, the new 
method is usable at least up to n « 30 (cf. ITable lj) . removing a significant 
obstacle for the development of practical random generators. 
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